Multiple defect model for non-monotonic structure relaxation 
in binary systems like Pd— Er alloys charged with hydrogen 
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In binary metallic systems like the Pd-Er alloys charged with hydrogen the observed structure 
evolution exhibits complex dynamics. It is characterized by non-monotonic time variations in an Er- 
rich fraction respect with an Er-poor fraction observed experimentally. The present paper proposes 
a qualitative model for this non-monotonic structure relaxation. We assume that the alloy have 
crystalline defects which trap (or emit) an additional amount of Er atoms. Hydrogen atoms into 
the alloy disturb the phase equilibrium as well as change the defects capacity with respect to Er 
atoms. Both of these factors lead to the spatial redistribution of Er atoms and cause the interface 
between the Er-rich and the Er-poor phase to move. The competition of diffusion fluxes in system 
is responsible for non-monotonic time variations, for example, in the relative volume of the enriched 
phase. We have found the conditions when the interface motion can change its direction several 
times during the system relaxation to a new equilibrium state. From our point of view this effect 
is the essence of the hydrogen induced non-monotonic relaxation observed in such systems. The 
numerical simulation confirms the basic assumptions. 
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I. INTRODUCTION 



For the last years we obtained a series of experimen- 
tal dataiiSiiii^ showing that palladium alloys (Pd-Mo, 
Pd-Ta, Pd-Er, etc.) mechanically strained and also 
charged with hydrogen exhibit non-monotonic time vari- 
ations in the phase structure during relaxation. Using 
X-ray diffraction techniques we studied alloys (Pd-Er 
in this paper) before and after the hydrogen charging. 
A palladium bulk was alloyed with Er of 8 atomic per 
cent then homogenized at temperature of 900 °C during 
24 hours and quenched. Further the specimen surface 
was mechanically polished, thereby the superficial layer 
of thickness of 2-5 /im became strained. Charging the 
specimen with hydrogen was performed electrolytically, 
in 4 % NaF solute twice distilled water under current flux 
of 80 mA/cm^ during 60 min. The X-ray diffraction pic- 
ture was monitored in real time by the automatic X-ray 
diffractometer using monochromatic CuKd irradiation. 
After charging with hydrogen the palladium specimen is 
kept under room conditions. 

At the first step we studied the Pd-Er alloy specimen 
after the polishing and before the chargingiiSiS,. AH the 
diffraction lines were of the Lorentz form in which the 
maximum location did not exactly correspond the face- 
centered cubic lattice but slightly differ to an individ- 
ual distance for each diffraction maximum. This points 
to the presence of anisotropic elastic strains in the X- 
ray reflecting superflcial layer that increase the spacing 
between the lattice planes parallel to the specimen sur- 
face. Also the diffraction maxima were slightly widened 
which indicate a large amount of dislocations in the re- 
flecting layer. An evidence of the two phase state was 
not fixed. Peculiarities of the diffraction picture we re- 
late to the bulk properties of the polycrystalline speci- 
men. The grain boundaries have a minor effect because 



of their small thickness. 

At the initial stage of relaxation directly after charging 
the homogeneous alloy decomposes on two phases and a 
large amount of hydrogen forms and restructs crystalline 
defects. With time hydrogen atoms degassing the alloy 
keeping the structure substantially non-equilibrium. A 
certain amount of hydrogen remains inside the system 
and bounds with Pd and Er atoms and defects that al- 
low spatial structure to evolve due to the interaction be- 
tween the nonequilibrium crystalline defects and the new 
phase clusters. Besides, a large amount of the crystalline 
defects accelerates diffusion of the solute atoms promot- 
ing the microstructure evolution. For example, the va- 
cancy concentration can attain 1-20 % inside metal- 
palladium alloys after charging with hydrogen (Pd-M- 
H systems )2i2424iS. That not only leads to the hydro- 
gen induced atom migration but also to forming metal- 
hydrogen- vacancy phase. 

In general the hydrogen charging changes the state of 
the Pd-Er alloy dramatically. Key events follow next. 

First, embedded hydrogen atoms increase the lattice 
parameter substantially. Under relaxation the elastic 
strain changes from the stretching to the contraction. 
Then it grows during the next two days and then goes 
down during the further eight days decreasing for 25 %. 
Hereafter, as time went up to one half year, no strain 
variations were detectedi*2i^. 

Second, the Lorentz form of diffraction lines now takes 
permanently the bimodal form. 

Third, under relaxation the profile of the diffraction 
line exhibits quasi-periodic variations. Two components 
of each diffraction maximum show oscillatory changes 
with respect to each other. That correspond to quasi- 
periodic variations in the relative volume, for example, 
of the phase enriched with Er atoms with respect to Er- 
poor phase. It should be noted that changing of the 
phase partition changes the Er concentration in it (for 
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FIG. 1: Experimental data for the Pd-Er alloy (8 at. % Er). 



example, the growth of Er-poor phase decreases Er con- 
centration in it) and, as a result, the difference between 
Er concentrations in two phases becomes greater. 

The entire phase structure of the Pd-Er alloy at differ- 
ent Er-concentrations is sufficiently complesii and con- 
tains many phases different in chemical composition. 
However we assume that corresponding phase transitions 
did not take place in the analyzed system because the 
mole fraction of Er atoms was sufficiently low, about 
10 %. The observed phase structure of the Pd-Er alloy 
evolves the solid solution of Er atoms in the Pd matrix 
and an Er-rich phase which should be related in struc- 



ture to the ErPdr intermetallic. Both of them may be 
modified with presence of H atoms. 

We explain observed phenomena like following. First, 
the initial elastic strains are due to defect -metal (D-M) 
complexes stretching the crystalline lattice by an image 
forces. Hydrogen change the strain sign by the conver- 
sion of the stretching D-M complexes into contracting 
H-D-M ones. Indeed, due to the high bond energy of 
the metal-hydrogen interaction H atom should be placed 
inside a defect and then pulls closer metal atoms sepa- 
rated by the defect. The defect effective volume becomes 
less than the volume of the matrix crystal cell. 

Second, the affinity of erbium for hydrogen is higher 
than that of palladium. Thereby, on one hand, the hy- 
drogen amplifies the Er-Er interaction. That decreases 
the stability of the Pd-Er solid solution and leads to 
decomposition of homogeneous state on Er-rich and Er- 
poor phases. On the other hand, lattice defects (disloca- 
tions, twin boundaries, etc.) contain in nearest volume 
a large amount of vacancies. In presence of hydrogen 
these regions should attract Er atoms. In other words, 
the lattice defects begin to play the role of Er sink. As 
a result the structure evolution of the Er-Pd alloy af- 
ter the hydrogen charging is governed by several factors 
which is responsible for the complex behavior of the sys- 
tem relaxation. Recently Albert Katsnelson et almi^ have 
proposed a macroscopic phenomenological model for the 
non-monotonic relaxation in the framework of synergetic. 
This paper presents a microscopic description of studied 
process. 



II. MODEL 

What is the essence of the observed phenomenon? We 
assume that homogeneous state of the Er-Pd alloy be- 
comes unstable during the hydrogen charging. A large 
amount of hydrogen atoms promotes decomposition on 
Er-rich and Er-poor phases and also changes the equilib- 
rium Er-concentrations of both phases. Also Era;H bonds 
are originated because of the higher affinity of erbium 
for hydrogen with respect to palladium. This process 
amplifies the effective Er-Er attraction and changes the 
phase instability threshold. That makes nonequilibrium 
the initially formed phase regions. The further relaxation 
is accompanied with the phase interface motion, leading 
to variations, for example, in the relative volume of Er- 
rich phase. Caused by Er-Er attraction redistribution 
of Er atoms takes place first in the close vicinity of the 
phase interface T. As a result, near the phase interface 
T the values Cq and cq of the Er concentration in Er-rich 
and Er-poor phases attain new equilibrium values c^{H) 
and co{H) depending certainly on the mean hydrogen 
concentration ch in the Er-Pd alloy at the time of the 
experiment. The values Cq{H) and cq{H) do not coincide 
with the Er concentrations inside the bulk just after the 
hydrogen charging. Thereby Er atoms diffuse away from 
the phase interface (or in opposite direction) causing the 
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FIG. 2: Schematic illustration of the Er distribution at the 
beginning of the relaxation process. 



interface T to move. With time the Er distribution near 
the the interface T becomes smoother and the speed of 
the system relaxation (interface velocity) becomes slower. 

The situation is complicated by crystalline defects 
(Fig. 121). On one hand vacancies and theirs complexes 
in Er-poor phase plays the role of Er traps. On the other 
hand defects inside Er-rich phase can be Er atom sources. 
Indeed, at the initial stage of the hydrogen charging de- 
fects can attract a large amount of Er atoms origining 
probably not only the pha se ErPdr but also the next 
phase ErPda (see Ref. [HI). In equilibrium the com- 
pound ErPdy is intermetallic with the narrow compo- 
sition interval inside which it exists. A large value of 
hydrogen affects the state of the compound ErPdy (or 
rigorously ErPd^Hj, where x w 7) enabling the phase ex- 
ist inside more wide composition interval. In this case 
ErPda phase can convert into the main Er-rich phase 
(ErPd^^Hj,) releasing extra Er atoms. We see originating 
of an Er source. All these defects, traps and sources of Er 
atoms, give rise to additional diffusion fluxes competing 
with the diffusion fluxes induced by the phase interface. 
Resulting flux affects the interface motion. Exactly this 
flux's competition is responsible for the non-monotonic 
relaxation. 

Let us discuss this effect in more details (with help of 
Fig. At the initial stage of the hydrogen charg- 

ing the chemical composition of the Er-poor phase (Er 
solid solution in Pd-matrix) consists of a higher value of 
Er concentration, cq, in comparison with the new equi- 
librium concentration cq{H) at finishing the hydrogen 
charging. This is due to H atoms in crystalline lattice 
amplify the effective Er-Er attraction. By the same rea- 
sons the opposite situation occurs in the Er-rich phase. 
In this phase the initial Er concentration is lower than 
the equilibrium concentration after the hydrogen charg- 
ing. Close by the phase interface T in the Er-poor phase 
the equilibrium values of Er concentration cq is constant 
during the motion of the interface T. It is due to the fast 



redistribution of Er atoms in area of the interface T. The 
resulting distribution of Er atom and the induced diffu- 
sion flux near the interface T is shown in Fig. |21 The 
velocity v of the interface obeys the expression 



(1) 



which depends only on diffusional coefficients D in each 
of phases and on gradients of Er-concentration close by 
the interface. In this formula the positive direction is 
shown by the axis of ordinates in Fig. [3 Ac is the drop 
in Er concentration c on two sides of the interface T, D 
and D* are the diffusion coefficients of Er atoms in Er- 
poor phase and Er-rich one, respectively, and, finally, the 
normal gradients V* c and V„c (taken at T) are directed 
diffusion fiuxes. We relate Er-rich phase to the ErPdj 
intermetallic which can exist only inside a narrow inter- 
val of chemical composition (V*c is small). Therefore 
we assume the diffusion flux inside Er-poor phase (Er 
solid solution in the Pd-matrix) dominates the interface 
motion. In other words, we approximate the rigorous 
expression ^ as 



w Ac w D\InC, 



(2) 



which focus our attention on Er-poor phase only. As 
follow from expression ^ the phase interface T should 
move inwards Er-poor phase (from left to right in Fig. 2). 
If this phase is homogeneous then the interface motion 
is monotonic (without changing direction). Its speed de- 
creases with time as the Er distribution becomes more 
and more uniform. 

The trap-defects in Er-poor phase change the situation 
dramatically. Evolution of the Er distribution in the de- 
pleted phase is illustrated in Figs. |3^a)-|3Id). Hydrogen 
activate some defects as Er atom traps. Near the defects 
Er concentration decreases and, as a result, takes a two- 
minima form shown in Fig.|2Ia). If capacity of the trap is 
sufficiently high then with time the Er distribution will 
take the form shown in Fig. Olb). In this case the Er 
concentration gradient V„c at the interface T changes 
the sign causing the interface T to move in the opposite 
direction. As time goes on the defect is filled with Er 
atoms, the Er distribution becomes again monotonic and 
the interface T moves in the initial direction. It should 
be noted that the source-defects in the Er-rich phase and 
the trap-defects in the Er-poor phase affect the interface 
motion in the same direction as it follows from Fig.|21 So 
we will consider only Er-poor phase. Er redistribution 
inside Er-rich phase cannot change essentially the results 
of our model. Its can change only the numerical results 
which are not principal for today. 

Mentioned above process shows changing direction 
of the interface motion two time, or shows one quasi- 
periodic fragment of the time variations in the phase par- 
tition. At the end of this fragment teh Er distribution 
becomes smoother (Fig. |2Id)) and the interface velocity 
is less than at beginning. The following quasi-periodic 
fragments of the non-monotonic relaxation is due to other 
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FIG. 3: Evolution of Er distribution inside Er-poor phase 
during the interface motion {ta < tt < tc < td)- 



defects in the Pd-matrix. We need to develop a concep- 
tion of their activation as traps or sources of Er atoms. 

Before the hydrogen charging defects expand Pd- 
matrix and the structure of the crystaUine defects should 
reflect this feature. Conversely after the hydrogen charg- 
ing defects squeeze the lattice. Under this conditions 
reconstruction of the defect should take a certain time 
before atoms rearrange themselves. It is quite reason- 
able that this process meets a certain potential barrier 
of the lattice elasticity nature. The elastic stress caused 
by other defects can make this barrier higher. But when 
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FIG. 4: Scheme of competing diffusion fluxes inside one do- 
main of the polycrystal. The process of non-monotonic relax- 
ation is kept by defects of different types activated at different 
time after hydrogen charging. 



nearest defects are reconstructed completely then the po- 
tential barrier is negated because of atoms arrange its 
structure. In this case reconstruction of the given de- 
fect becomes more easy. From this point of view we as- 
sume that the defect will be activated as the trap of Er 
atom only when the preceding defect is practically filled 
with Er atoms. In other words, the defect activation 
during the interface motion proceeds like the domino ef- 
fect. Each act of the defect activation and filling with Er 
atoms corresponds to one quasi-periodic fragment of the 
non-monotonic relaxation. Besides there are a number 
of crystalline defects of different size and kindiS which 
can make the relaxation even more complex. The larger 
is a defect, the higher is the potential barrier of atom 
rearrangement and, so, the more is the time needed for 
its reconstruction. Thereby small defects should be acti- 
vated first and only then large defects can come into play 
(see Fig.01). The more is the distance between the defect 
and the interface (or the bigger is the defect), the more is 
the time of Er redistribution near the interface. There- 
fore at the further stage of the relaxation the duration of 
quasi-periodic fragments is expanded remarkable. 

In the present paper we will consider only characteris- 
tics of defects which are needed us for qualitative analy- 
sis. Our aim is point reader to key moments of structure 
evolution in Pd-Er-H system. Quantitative description 
is the work of future. Now we can formulate the gov- 
erning equation simulating the interface motion and the 
non-monotonic relaxation. 
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III. GOVERNING EQUATIONS 

Taking into account aforesaid we will consider the dif- 
fusion of Er atoms inside Er-poor phase (regarded as the 
half-axis z > zq-) bounded by the phase interface T. We 
write the diffusion equation for the atomic concentration 
c of Er atoms in this phase as 



dc 
dt 



D 



dz"^ 



I 



^ qiJi{qi^i)S{z - Zi) . (3) 



i=l 



Here D is the diffusion coefficient of Er atoms (which 
is constant in this article), the trap-defects are approx- 
imated as the (5-like sinks placed one from other at the 
equidistant distance {zi = iz^^}°Zi, ^^'^ the prime means 
that the sum runs over all the defects located in the Er- 
poor phase, Zi > zq-- A 5-sink is characterized by the real 
physical size I of a trap-defect, by the average time Ttr 
of trapping an Er atomiSi, and by the defect capacity q, 
i.e. the number of free seats for Er atoms at the current 
moment of time. 

At the initial time t ^ the capacities of all the de- 
fects have the same value go- As time goes on the active 
defect i traps Er atoms decreasing the capacity qi. We 
assume defect i to be activated when the capacity of the 
preceding defect i — 1 drops down the threshold Qc = OcQo 
where 0c ^ 1- This behaviour is described by the acti- 
vation function for i > 2 



Ji{qi-i) 



0, if > (jc, 

1, if qi^i < qc- 



(4) 



For the first defect, by definition, Ji ~ 1, because we 
assume it becomes active by the phase interface directly. 
The defect capacity obeys the equation 



^ = — c{zi)qiJi{q.i^i) 

dt aTtr 



(5) 



where a is the spacing of the Pd-lattice. 

The phase interface T initially is located at zero point 



ZT\t=0 



0. 



(6) 



Its further motion is governed by the diffusion flux of Er 
atoms inwards the Er-poor phase and obeys the equation 
(see formula ^) 



dof dzr D dc 



dt 



Ac dz 



(7) 



where the value Ac is shown in Fig.[2| In writing expres- 
sion (0 we ignored the effect of diffusion flux inside the 
Er-rich phase on the interface motion. Besides when the 
phase interface crosses the point Zi (so defect i appears 
in the Er-rich phase) we assume it not affect the interface 
motion further. This is justified by the fact that the in- 
terface T can pass through the defect location if only it is 
filled with Er atoms and, thus, invisible to the interface. 



At the interface T the Er-concentration takes the new 
quasi-equilibrium value which formed by fast local Er 
redistribution near T (Fig.|21): 



c{zt) = Co - AcH 



(8) 



where cq is the initial Er-concentration inside the Pd- 
matrix, i.e. 



c{z,t)\t=Q ^ Cq for z > 0. 



(9) 



The system of equations (O, ©, with the boundary 
condition (|SJ| and the initial conditions ©, forms 
our model of the phase partition evolution or, generally 
speaking, the non-monotonic relaxation. 



A. One quasi-periodic fragment of non-monotonic 
releLxation 

Let us, first, find the main features of the model and 
write parameters (in dimentionless form) governing the 
structure relaxation. To do this we consider the begin- 
ning of the interface motion when the only one nearest 
defect is active. In the next section we will study the 
array of defects. 

For convenience we attach our frame of reference to the 
moving interface. We will use in formula dimensionless 
time T — {Dt) / (z^j.)'^ , coordinate ^ = (z — zr)/z^j., nor- 
malized Er concentration rj — c/cq, the defect capacity 
^ = 91/90 and the dimensionless time-dependent coordi- 
nate ^tr(T) = (z°j. — zr{t))/z^j. of the first defect in the 
moving frame. In these terms the aforestated system is 
converted to the following equations: 



Dt 



(11) 



Equations (|10|l . must be completed by the boundary 
conditions for r > 0: 



Ac dftr(r) 

Co 



dr 



V{^,r)\^=0 = 1- 
and the initial conditions: 

^KC,^)l.=o = i^ ^(t)|..o = i^ 

for ^ > 0. Following quantities: 

Ac AcH ^ _ qplz^^ 
Co ' Co ' Dtu- ' 



AcH 
Co 



(12) 
(13) 



6r(T)k=0-l (14) 



A = 



„ ,0 
go a 



(15) 



form the dimensionless parameters governing the model. 

f2 is one of the key parameters of the model. If 51 ^ 1 
(which corresponds to small capacity of the defect go or 
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to very fast diffusion) then the trap-defect not affect the 
interface motion. Lets to be of order unity. At the 
initial stage when t ^ 1 (i.e. t ^ (z°j.)^/_D) the only Er 
redistribution near the interface T governs its dynamics. 
In this case the latter term on the right-hand side of 
Eq. (|10(l can be ignored and the system of equation pO|l . 
(|12|I - H14() has a self-similar solution: 



Co 



1 - erf 



where 'i?o is the root of the transcendental equation: 



(16) 



/^^oe"^" [1-erf (z9o)] 



AcH 



(17) 



The corresponding value of the dimensionless interface 
velocity -i^ = vz^^/D is written as i) — ^q/^/t. For nu- 
merical experiment we take next values: Ac/co = 1.5 
and Ach/co = 0.5. In this case the value of do is ap- 
proximately equal to -do ~ 0.2. For the given values the 
interface should reach the first defect (with ^ 1) in 
time T « 10. 

Now let us activate the defect with ft = 1. As follows 
from Eq. (|10|1 for a sufhciently short time t ~ all the 
Er atoms near the trap-defect will be trapped. Therefore 
in the vicinity of the defect Er distribution takes the form 
shown in Fig. Ola). The general form of Er distribution 
near the defect remains unchanged until either the inter- 
face T comes close to the defect or the defect is filled with 
Er atoms. Certainly, the distribution width will increase 
with time as ^/T. From Eq. for fl^'^ <C t <C 1 Er 
distribution near the defect can be approximated as 



r]iC,T) = erf 



1$ ~ gtr 



(18) 



Within the same time interval the dimensionless rate of 
trapping Er atoms by the defect is estimated as 



dr] 



(19) 



Because of diffusion dithering at time r 1 Er distribu- 
tion will get the form shown in Fig. I^^b). Such a distri- 
bution has a negative value of the boundary gradient, so 
the interface T moves in the opposite direction (accord- 
ing to expression H12|l or Q). As time goes on further the 
defect will be filled with Er atoms. The Er concentration 
in area of the defect will grow and the defect will affect 
weaker the interface motion. The corresponding form 
of Er distribution is illustrated in Fig. EJc). As follows 
from Eq. Ijlll) and estimate (|19|l the time of the defect 
filling is about r ~ A~^. After that the defect cannot 
affect the interface motion. Er distribution tends again 
to the self-similar form (|16|l . For observer the interface 
velocity changes the sing for the second time, and the 
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FIG. 5: The interface dynamics described by the system l|10|l - 
11411 for different values of the defect capacity. The fragment 
(a) shows the dimensionless interface velocity vs time and 
the fragment (b) illustrates the time dependence of the inter- 
face displacement. The abscissa exhibits time at root-square- 
law scale. In numerical simulation we used: Ac/cq = 1.5, 
Ach/co = 0.5, and fl = 1. The values of the parameter A are 
pointed out at the curves. 



interface motion returns to the initial direction until the 
next trap-defect is activated. 

The aforementioned scenario describes one quasi- 
periodic fragment of the non-monotonic structure re- 
laxation. For numerical analysis of (|10|I - H14(I we take 
Ac/co = 1.5, Ach/cq = 0.5, fl = 1 and then vary 
parameter A (in non-mathematical words we study de- 
fects of different capacities). The results are shown in 
Fig. [SI When the capacity of the defect is sufficiently 
high (A < 1) wc obtain an essential backward interface 
motion. Otherwise, the trap-defect will be filled before 
the induced diffusion flux could affect the interface mo- 
tion. 

When the first defect is practically filled with Er-atoms 
the next defect will be activated. That should give 
rise to a similar one quasi-periodic fragment of the non- 
monotonic relaxation, which is the subject of the next 
section. 
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B. Multi defect effect on the interface dynamics 

In this section we will use the same dimensionless time 
r, coordinate ^, the normalized Er concentration 77 and 
defect capacities {Oi = qi/qo}- In these terms the initial 
system of governing equations is reduced, first, to the 
diffusion equation written in the form 



where according to Q for i > 2 



(20) 



0, if > 

1, if e,^i < 



and Ji = 1. 



In expression H20|l the value d = vz^^/D^ as before, is 
the dimensionless velocity of the interface T and the 
sum runs over all the defects located inside the Er-poor 
phase (whose coordinates ^^(t) > are positive in the 
frame attached to the moving interface T). By definition 
^^(t) = i — zr{T)/ z^^ all the quantities ^^(t) changes with 
time at the same rate 



dUr) 
dr 



(21) 



Second, the defects filling is described by the system of 
equations 



^ = -A17,?(e.(r),r)0,J,(0,_i) 
dr 



(22) 



Third, the interface velocity d relates with the boundary 
condition as ((T^ : 



9r,(^,T) 



Ac 



Co 



(23) 



The boundary condition (|13|) holds. The initial condi- 
tion H14() should be replaced by the following one 

^(C,r)|,^o = l, 0.(r)|,^o-l, ^,(r)|.=o = * (24) 

for ^ > 0. 

The system of equations has been analyzed numeri- 
cally for the same values of the parameters Ac/cq = 1.5, 
Ach/cq — 0.5 but with fixed the value of the parame- 
ter A equal to A = 0.5. In other words here we fix the 
defect capacity and analyze the interface dynamics de- 
pending on Er diffusivity and the threshold value 9c. At 
first we set Oc = 10"^ and Vl — 2. The corresponding 
dynamics of the interface T is shown in Fig. |3 Similar 
non-monotonic structure evolution is observed for O ~ 1. 
When the parameter takes sufficiently large values (i.e. 
if the Er diffusivity D is small enough) so that ^ 1 
then the trap-defects affect the interface motion as one 
trap but of much higher capacity. For example. Fig. [7| 
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FIG. 6: Dynamics of the interface T described by the sys- 
tem I|2()|l - lj24|l in which the interface motion is affected by 
four trap-defects activated sequentially. The fragment (a) 
shows the dimensionless interface velocity vs time and the 
fragment (&) illustrates the corresponding time dependence 
of the interface displacement. The abscissa exhibits time at 
root-square-law scale. In numerical simulation we used the 
following values of parameters Ac/co = 1.5, Ach/cq = 0.5, 
A = 0.5, = 10"*, and set Q. = 2. 



exhibits the interface dynamics for U, — 10. In the given 
case the interface motion looks like one quasi-periodic 
fragment but prolonged substantially. For 17 ^ 1, i.e. 
when the Er diffusion is sufficiently fast the defects will 
be filled also fast and they have no considerable effect on 
the interface motion. 

The specific details of the interface dynamics depend 
on the threshold Oc of the defect activation. In particular, 
for Oc = 10~^ we observe the interface dynamics similar 
to Fig. for Vl — 1 whereas the value = 2 got us to the 
interface motion of the same form as shown in Fig. 

It should be noted that in Fig.[5l|7|the abscissa is a time 
at root-square-law scale in order to emphasize the oscil- 
latory behaviour of the non-monotonic structure relax- 
ation. In reality each following quasi-periodic fragment 
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FIG. 7: Dynamics of the interface when the Er diffusion is 
fast. It show the interface displacement vs time predicted by 
the system 1201 - 12411 when the interface motion is affected 
by four trap-defects activated sequentiaUy. The abscissa ex- 
hibits time at root-square-law scale. In numerical simula- 
tion we used the following values of parameters Ac/co = 1.5, 
Ach/co = 0.5, A = 0.5, 0c = 10~*, and set Q = 10. 

is more prolonged than the preceding one. A similar be- 
haviour was found experimentally (see Fig.^). However, 
the experimental curve exhibits quasi-periodic time vari- 
ations at logarithmic scale rather than root-square scale 
for the theoretical curve. We explain this discrepancy by 
taking into account different size of defects. When the 
larger a defect is, the later it comes into play and the 
longer time it will be filled with Er. 

IV. CONCLUSION 

We have proposed a mechanism which responsible 
for the nonmonotonic structure relaxation observed in 



metallic binary systems like the Pd~Er alloys after hydro- 
gen charging. The observed quasi-periodic alterations of 
the diffraction line profile are related with time variations 
of volume and concentration of two phase (Er-rich and 
Er-poor) arose in the alloy after hydrogen charging. We 
developed the model in which we consider active crys- 
talline defects that trap additional amount of Er atoms. 
The key point of the model is disturbance by hydrogen 
the phase equilibrium as well as change the defect capac- 
ity with respect to Er atoms. Both of these factors lead 
to the spatial redistribution of Er atoms causing the in- 
terface to move between the two phases. Because of the 
induced diffusion fluxes these factors individually would 
move the interface in the opposite directions. As a result 
their competition is responsible for non-monotonic time 
variations, for example, in the relative volume of Er-rich 
phase observed experimentally. 

We have found the conditions when the interface mo- 
tion can change its direction several times during the 
system relaxation to a new equilibrium state. The latter 
effect is, from our point of view, the essence of the hydro- 
gen induced non-monotonic relaxation observed in such 
systems. The numerical simulation confirms the basic 
assumptions. 

It should be also noted that the present paper pre- 
tends only to a qualitative description of the observed 
non-monotonic structure relaxation. Actually we have 
singled out the basic feature of such systems that is re- 
sponsible for this complex process. To compare the ex- 
perimental data with results of numerical simulation a 
more sophisticated model should be developed. 
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